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Abstract 

We generalize the KPZ equation to an 0(3) N = 2j + 1 component model. In the 
limit N — > oo we show that the mode coupling equations become exact. Solving these 
approximately we find that the dynamic exponent z increases from 3/2 for d = 1 to 2 
at the dimension d ~ 3.6. For d = 1 it can be shown analytically that z = 3/2 for all j. 
The case j = 2 for d = 2 is investigated by numerical integration of the KPZ equation. 
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Many growth models have been studied in recent years both analytically and numerically 
A widely used description of the broad scale features of such growing surfaces is a non-linear 
Langevin type equation proposed by Kardar,Parisi and Zhang (KPZ) U : 



J^(r, t) = z/V 2 + A(V0) 2 + V (r, t). (1) 



It appears to describe the surface profiles generated in the Eden model ||] and by ballistic 
deposition [[| . It can also be easily mapped onto the directed polymer problem in a random 
potential 0] . It describes the time evolution of a single valued height parameter 0(r, t) (ie 
with no overhangs or voids) for a growth process onad dimensional substrate. The equation 
reflects the competition on mesoscopic length scales between surface tension smoothing forces 
z/V 2 0, the tendency for the growth to occur preferentially in the direction of the local normal 
to the surface, represented by the term A(V0) 2 , and the noise which is Gaussian, such that 

(r ] (v,t)r ] ( r ',t')) = 2D6 d (r-r')6(t-t'). (2) 

The objective is to characterize the form of the surface. A step towards this is to evaluate 
the two-point correlation function, 



C(k u) - Mk,^'(kW)) (3) 
° lk ' UJ) ~ (27r) d +W(k - k')6(u - u>) {A) 



and the response function in frequency space, 



fYk \ l / <90(k, u) x 



There are two distinct regimes in the scaling limit for d > 2. In the weak coupling limit, i.e. 
for A < A c , the behavior is governed by the A = fixed point. In the strong coupling regime, 
A > A c , due to the non-linearity the correlation and response functions take on the non-trivial 
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scaling forms, 



Gfr") = (5) 



where x is related to z via the scaling relation \ + z — 2- In the weak coupling regime z = 2 
and is independent of d. For d < 2 only the strong coupling regime exists |5|, § . Knowledge of 
the exponent 2 as a function of the dimensionality d, and the scaling functions, would mean 
the system was essentially understood. 

Some workers 0, |8|, H believe that above some finite critical dimension d c , z = 2 for both 



the strong and weak coupling regimes. This is in opposition to numerical work |U| which 
suggests that the upper critical dimension is infinite, although crossover effects might mask 
the true value of z for the small systems studied in higher dimensions. 

A major theoretical difficulty is that in the strong coupling regime the perturbation series 
in A about A = cannot be summed self-consistently in terms of just response and correlation 
functions because of vertex correction graphs which renormalize the non-linearity as shown 
in Fig. 1. (The perturbation formalism we are using is that of Ref. ||.) However as at each 
order in A graphs containing vertex corrections scale in the same manner as graphs without 
vertex corrections, it is tempting to ignore all the vertex correction graphs completely and 
sum the graphical series Fig. 2 self-consistently , leading to the mode coupling equations ||11|| . 
These equations are expected to give qualitatively correct values for the exponent z yet they 
remain an uncontrolled approximation. Some modification of the KPZ equation for which 
mode coupling was exact would place the approximation on a surer footing. 

A widely used method in the study of critical phenomena is to increase the number of 
components of the field to N, where it is often found that the new model is exactly solvable 
in the infinite component limit and that a systematic 1/N expansion may be developed. It 
is possible to perform such generalizations in many ways. Kraichnan and Chen fl3| in their 



formulation of the DIA equations for turbulence took a model where the vertex factors had a 
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random amplitude ±1/JV. As in the turbulence problem considered by Mou and Weichman 
14| , we use the approach of Amit and Roginsky [[15j and generalize the KPZ equation to an 



iV = 2j + 1 component field which forms the basis for an irreducible representation of the 
0(3) symmetry group, 

|> m (r, t) = z/V 2 m + 4=J2 ^ n V0;V0 n + Vm (r, t). (6) 

The non-linear term in the KPZ equation represents the coupling of two fields to produce a 
field of index m, m — —j, — j + 1, . . . , +j. The noise has the form, 

( Vm (r,t) V * m ,(r',t')) = 2D5 m>m ,5 d (v - r')5(t - t') (7) 

with rf m = (-l) m r]_ m . 

For the form of the generalized KPZ equation to be independent of the representation it 
must be unchanged by the transformation (fi m — > <p m i = 7£™/(ii)(/> m given u G 0(3) , 1Z™, being 
NxN matrices which form an irreducible unitary representation of 0(3). The coupling tensor 
A 1 ^ 1 must satisfy 

Kf' = KM") K'H ftf («) ■ (8) 

The tensor A l £ is therefore the Clebsch-Gordan coefficient (j, n;j, l\j, m). The angular mo- 
mentum j must be an even integer for non-zero A 1 ^ . The case j = (N = 1) is the scalar 
KPZ equation fll]). 

We shall now consider the large j limit. The structure of the graphs is the same as that of 
the 3 field theory studied by Amit and Roginsky [15]. Taking over their results regarding the 



asymptotic form of Clebsch-Gordan coefficients it can be seen that graphs with one or more 



vertex corrections have magnitude 1/N or less [16] relative to other graphs with no vertex 
corrections of the same order in A. Thus all the vertex corrections are negligible as iV — > oo 
and the graphs for the correlation and response functions may be summed self consistently 
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by performing the summations over the component labels using the orthogonality relation , 



Y.UJ;j,n\j,m)(j,l;j,n\j,m') = 5 m ,m<- (9) 
n,l 

The correlation function C(k,u) = (0 m (k, c<j)0* (k, uj)) and the response function related 
to (d(f> m (k, uj)/dr] n (k, u>)) vanish unless n = m , and for n = m they are m independent as a 
consequence of the overall rotational invariance of the equations. 

The mode coupling (MC) fll| equations can be written as, 



C{k,u) = 2D\G{k,u)\ 2 



{27l) d+1 

G _1 (k,w) = G \k,uj 
4A 2 



/+oo f+OO 
rf/i / d d q (k + .k_) 2 C(k + ,w + )C(k_,cj_) (10) 
-oo J —oo 



(27T) 



d+1 



/+oo f-\-oo 
dfi d d q (k + .k_)(k + .k)C(k + ,cj + )G(k_,^_) (11) 
-oo J-oo 



where k± = y ± q and u;± = y ± /i. The first term in Eq. (|T0| ) can be dropped in the 
scaling limit provided z < (d + 4)/3. We find this condition is satisfied by our solution. 
In Eq. (|Tl"|), Go 1 (k, uj) = uk 2 — iu, but only the term — iuj is relevant in the scaling limit. 
Inserting the scaling forms for the response and correlation function into the mode coupling 
( p~Qj . p~T| ) equations yields the scaling relation x + z — 2. To calculate the exponent z from 
the MC equations we need to input the exact scaling forms for the response and correlation 
functions in the right hand side of Eq. (|i~0"yi~lD where consistency will only be achieved for the 



exact value of z. A previous attempt by Bouchaud and Cates [|17] to solve the MC equations 
(approximately) for the exponent z by this consistency requirement took the response function 
in k, t space to be a simple exponential decay , 

G{k,t) =exp{-k z t) 6{t) (12) 
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and the correlation function to be, 



C(k, u) = Vk 6z - d -*\G{\ii (13) 

where T> is a constant. We shall use instead as the starting point of our calculation the work 
of Hwa and Frey |ISJ in d = 1 on the form of the scaling function. By solving the mode 
coupling equations numerically they showed that 

G(k, t) « exp{-ak 2 t 2/z ) 6{t) (14) 



was a good approximation for the response function [p0[| . The parameter a is dependent on 
the units used and may be scaled out by the transformation k — > k/-y/a. Since in d = 1 and 
only in this dimension the correlation and response functions satisfy a fluctuation dissipation 
theorem (FDT)0, 

C(k, W ) = l(G(k,w)+G*(k,u;)) J (15) 
we take the form of the correlation function to be 



^(M = F ^exp(-k 2 |tr) (16) 



where B is some arbitrary constant that depends on A. Note that the correlation function is 
an even function of uj and that it satisfies the FDT for d = 1 and z = 3/2. 

We shall assume that the form of the scaling functions do not vary strongly with the 
the substrate dimension d. We work away from d = 1 using the same choice for the scaling 
functions. The Gaussian form of the response and correlation functions allows most of the 
integrations to be done analytically. Matching both sides of the fourier transformed versions 
of (PT|) and (11) at u = allows z to be determined ||22|| . This matching procedure is an 



arbitrary choice out of many schemes available but it is hoped that the exponent z does not 
depend significantly on the details of the matching. 

6 



The form of the variation of z with d we find to be qualitatively similar to that of other 
self-consistent treatments Jl7|, jlSf . We obtain the exact result z = 3/2 for d = 1. Our estimate 
for z is 4/3 as d — > 0. (This is the exact value for z for any j as d — > 0.) The numerical 
value of the exponent was 2 = 1.662 for d = 2 which is close to the value found by Bouchaud 
and Cates [17], 2 = 1.67, but still some way off the value obtained from simulations for j = 
23 where z ~ 1.614 . The exponent z reaches its weak coupling value z = 2 at d ~ 3.6 



which is similar to the values found in Ref. |17j d c ~ 3.75 and significantly larger than the 
value of Schwartz and Edwards [|18| d c « 3.25. It should be noted that the form of the our 
scaling functions is identical to that of Bouchaud and Cates for z = 2 so d c , the upper critical 
dimension, should be the same in both approaches. We do not understand this discrepancy. 

The variation of z with d thus seems to be only weakly dependent on the assumed form 
of the scaling functions. The crucial question remains whether this apparent upper critical 
dimension is an artifact of using inappropriate approximations for the scaling functions. We 
hope in future work to develop a systematic procedure of matching the derivatives of the 
mode coupling equations at uj = . This will provide better approximations for the form 
of the scaling functions and more accurate estimates of z away from d = 1. It may also be 
possible to solve for the scaling function numerically (as Hwa and Frey did in d = 1) for 
general dimension d. 

The strong coupling value of the exponent z might be expected to vary continuously 
between j = (the scalar KPZ) and the j — > 00 limit (in which case mode coupling is exact). 
We have investigated the case of j = 2 which is an intermediate value. For j = 2 the field <j) m 
has 5 independent components. It is analogous to the case of a traceless symmetric 3x3 tensor 
order parameter. We calculated a rotationally invariant characteristic width of the interface 
at equal times, 

/ 24 m=2 \l 

W{t)= I— ( £ 2 m (r,t)CM)>j (17) 
by numerically integrating the KPZ equation in d = 1 and d = 2 using a discrete grid method 
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\24\. The width calculated at equal times for a finite size system of width L scales as, 



by analogy with the scalar KPZ equation pSfl . For short times W(L,t) ~ t x ^ z . We evaluated 
x/z from the slope of a log-log plot of W(L,t) versus t calculated at each time step 5t for a 
system with for d = 1, L = 20, 000, g = 10 (g being the magnitude of the effective coupling 
g = 2\ 2 D/v 3 ) and 5t = 0.001 to obtain, \jz = 0.33 ± 0.02. This value was found to be very 
robust for different g and St and is in agreement with the scalar KPZ equation in d = 1. 

The exponent x is determined by the usual procedure of starting from a flat surface of 
linear size L and growing it until the fluctuations in the height are saturated by finite size 
effects where W(L) ~ L x . W(L) was measured for system sizes L = 5, 10,20,40,80, 160 in 
d — 1. We plotted \og(W (2L) 2 — W(L) 2 ) versus log(L) to eliminate any constant correction 
to scaling. By least squares fitting of the data we obtained x — 0.50 ± 0.01. The value is the 
same as that of x f° r the scalar KPZ equation. 

We attempted to calculate the exponents by the same method for the case d = 2. The 
calculation was plagued by an instability where the width W(L,t) diverges after a certain 
period of time tj. This time was found to be only weakly dependent on the time discretization 
used. Smaller time steps 5t only delayed the onset of the instability by a small amount. 
Furthermore ti was found to decrease as A was increased. For times less than ti a region 
exists in which the \ogW(L,t) versus logt plot has a constant slope but this slope was found 
to rise as A was increased. This may be due to some crossover effect from weak to strong 
coupling. The maximum slopes are consistent with the exponent ratio x/ z ~ 0.16 — 0.2 giving 
z fa 1.67 — 1.72 . These difficulties are the same as those found for the scalar KPZ equation 



in d — 3 |26| , but the instability is more pronounced for j = 2 making reliable calculation 
of the exponents difficult. The instability may be due to the discretization of time and space 
used in the numerical integration, but further work is needed to substantiate this idea. 

Using the methods of Huse, Henley and Fisher H we have found the stationary point so- 
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lution in d = 1 of the equivalent Fokker-Planck equation for the N component KPZ equation. 
From this we see that a fluctuation-dissipation theorem exists in d — 1 for all even j, and that 
z = 3/2. To observe the effect of varying j on the exponent z the system has to be studied in 
d = 2. As we have seen, this presents a formidable numerical challenge. We speculate though 
that the generalized KPZ equation for j = 2 might actually describe a situation of physical 
interest e.g. the time evolution of the orientational order of a surface grown by aggregation 
of ellipsoids. 
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Figure Captions 

Fig. 1: An example of graphs which renormalize the value of the non-linear vertex term in 
the perturbation series for the correlation function. A. In the large N limit the contribution 
of these graphs is down by a factor of 1/N. The circle denotes the noise term 2D. 
Fig. 2: The graphical series neglecting vertex correction graphs for the correlation C(k, u) 
and response G(k,u) respectively, which give the mode coupling equations. 
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